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1. Introduction 



String-like configurations of classical background fields can polarize the vacuum of 
quantized matter fields, resulting in the emergence of a nonzero vacuum expectation 
value of the energy-momentum tensor in spatial regions where the classical field is 
zero. A particular example of a string-like configuration is a tube of the magnetic 
flux lines, that is formed inside a current-carrying long and thin solenoid. But a more 
general string-like structure is the Abrikosov-Nielsen-Olesen vortex (TJ |2] arising as a 
topological defect in the aftermath of a phase transition with spontaneous breakdown 
of a continuous symmetry; the condition of its arising is that the first homotopy group 
of the group space of the broken symmetry group be nontrivial. The vortex, apart from 
its flux, is characterized by nonzero energy distributed along its axis, which in its turn, 
according to general relativity, is a source of gravity. It can be shown [3j HI [5] that 
this source makes space outside the vortex to be conical with the deficit angle related 



Induced vacuum energy-momentum tensor 



2 



to the energy density per unit length of the vortex. Since the squared Planck length 
enters as a factor before the stress-energy tensor in the Einstein-Hilbert equation, the 
deviations from the Minkowskian metric are of the order of squared quotient of the 
Planck length to the correlation length, the latter characterizing the thickness of the 
vortex. For superconductors this quotient is vanishingly small and effects of conicity are 
surely negligible [Bj; hence, the appropriate vortices known as the Abrikosov ones are 
much the same as the solenoidal magnetic vortices (with one important distinction that 
the flux of the former is quantized). However, topological defects of the vortex type arise 
in cosmology under the name of cosmic strings [7J [8] . Cosmic strings with the thickness 
of the order of the Planck length are definitely ruled out by astrophysical observations, 
and there remains a room for cosmic strings with the thickness which is more than 3.5 
orders larger than the Planck length, that corresponds to the deficit angles bounded 
from above by the value of 2.5 x radians (see [9]). Cosmic strings associated with 
spontaneous breakdown of global symmetries are called global (or axionic) strings; they 
are characterized by zero flux. Cosmic strings associated with breakdown of gauge 
symmetries are called gauge (or local) strings; they are characterized by nonzero flux. 

Cosmic strings can have various astrophysical effects, in particular, they serve 
as plausible sources of detectable gravitational waves (TOj [TTJ [12], gamma-ray bursts 
[13], high-energy cosmic rays [H]. Although the primary role of cosmic strings in 
the formation of the large-scale structure of the Universe is ruled out by observation 
of cosmic microwave background by COBE and WMAP [151 H5]> cosmic strings can 
produce a mechanism for the generation of the primordial magnetic field in the early 
Universe [TTJ HH [19]. The interest to cosmic strings is augmented in the last decade 
owing to the finding that they emerge at the end of inflation in the framework of 
superstring models [201 EE] (see also reviews [22J [231 12H 125]). The details of the cosmic 
string formation mechanism can be quite different, either fundamental (super) strings or 
Dirichlet branes are stretched to a macroscopic size, but such structures arize almost 
inevitably in supergravity models with large extra dimensions [26]. In this respect it is 
appropriate to consider a generalization of a cosmic string in 4-dimensional space-time 
to a codimension-2 brane in higher- dimensional space-time. Such a brane appears also in 
the construction of brane- world cosmologies with two extra dimensions (see [271 ES [29] 
and references therein), as well as in the study of quantum aspects of the black hole 
physics, playing a key role in the formation of chains of Kaluza-Klein bubbles with 
black holes [301 13T] . Therefore, the properties of the vacuum of quantized matter in the 
background of a codimension-2 brane deserve investigation. 

A study of the vacuum polarization effects in the cosmic string background has 
a long history. The induced vacuum energy-momentum tensor was considered for the 
cases of global [321 [331 EI] and gauge [351 EH] strings; in the above the quantized fields 
were assumed to be massless, whereas the results for the nonvanishing mass were either 
incomplete or of restricted use, see [371 EH EH]. The effects of nonvanishing mass were 
exhaustively studied for the case of a negligible deficit angle only [10] (see also [El 02]). 
In this respect it should be noted that the induced vacuum current in the cosmic string 
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background is known both in the massive [TS] and massless [13] cases. 

The aim of the present paper is to obtain the energy-momentum tensor which 
is induced in the vacuum of the quantized massive scalar field in the cosmic string 
background. In the next section, the vacuum expectation value of the energy-momentum 
tensor is in general defined with the use of the point-splitting regularization. Green's 
function in the cosmic string background is obtained in section 3. The induced vacuum 
energy-momentum tensor is obtained in section 4. Some limiting cases of our result are 
discussed in section 5, while the case of the massless quantized scalar field is considered in 
section 6. The results are summarized in section 7. We relegate the details of derivation 
of Green's function in the cosmic string background to the appendix. 

2. Energy-momentum tensor: point-splitting regularization 

The energy- momentum tensor for the quantized charged massive scalar field ^>(x) is 
given by expression [HJ US] (for a review see [46J ) 

T» v {x) = l - [Vtt + (aO, V*(z)] + + \ [V^ f (x), Vtf (x)] + 

+~<T [(□ + m 2 + £R)*\x), *(*)] + + i<T [*t( x ), (□ + m 2 + + 

+(£- l)g*TL [&{z), *(x)] + -e(V t V + i2^) [&(x), 9(x)] + , (1) 

where V M is the covariant derivative involving both affine and bundle connections, 
n = V M V M is the covariant d'Alembertian operator, g^ u is the space-time metric tensor, 
R^ v is the Ricci tensor and R = g^R^ is the scalar curvature of space-time, £ is the 
coupling constant of the scalar field to the space-time scalar curvature. The product of 
the field operators at the same point (x) is ill-defined, and this leads to the divergence 
of the formal expression for the vacuum expectation value of the energy-momentum 
tensor, (0\T^ u (x)\0) . To regularize this divergence, one uses the point splitting for the 
operator product in (1): 

T^(x; x') = D» v (x; x') [T*(x)tf(x') + T*(x')&(x)] , (2) 

where T is the symbol of chronological ordering in the operator product and 

D^(x; x')= Q-^ (V^ + V^,) 

+&T + m 2 + £R X ) + ig^{U xl +m 2 + £R X ,) 

+2 U - ~) sT (gp*vm -m 2 - \iR x - l -ZR x ^j 

The vacuum expectation value of tensor T^ v (2) is well defined. Further, if field 
operator \1/ obeys the Klein-Gordon equation, then the vacuum expectation value of 
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the chronological operator product in the right-hand side of (2) is related to Green's 
function of the Klein-Gordon operator, 

(0\TV(x)tf(x')\0) = -iG(x; x'), (4) 

where 

(D x + m 2 + £R X ) G{x; x') = (n x , + m 2 + £R X ,) G(x; x') = —5{x-x'){h) 

V 9 

and g = |detg^2,|. Green's function is decomposed as 

G(x; x') = G(x; x') {S ) + G(x; x') {R) , (6) 

where the first term is singular, whereas the second one is regular in the coincidence 
limit, x' — > x; the divergence of the vacuum expectation value of tensor T^ u (1) is due 
to the contribution of G(x; x')(s) + G(x'; a?)(s) to tensor T^ u (2). Thus, the physically 
meaningful (renormalized) vacuum expectation value of the energy-momentum tensor 
is obtained by taking the coincidence limit after the subtraction of the dangerous 
(divergent in this limit) part: 

t^(x) = -i Urn D^(x; x') [G{x; x') (R) + G{x'; x) (R) ] . (7) 

This is a quite general scheme which is only necessary but not sufficient to define 
unambiguously the renormalized value. As we shall see in the next section, to get the 
renormalized value in the case of the cosmic string background, it suffices to subtract 
the appropriate value in the case of the trivial (Minkowski) background. 



3. Cosmic string background: Green's function 

As is already mentioned, a cosmic string associated with spontaneous breakdown of a 
gauge symmetry is characterized by a gauge field with the strength which is directed 
along the string and is non vanishing inside its core; the total flux of the strength is 
denoted in the following by $. Space-time outside the string core is the conical space- 
time where a surface orthogonal to the string axis is isometric to the surface of a cone; 
the deficit angle is 8irGfi, where G is the gravitational constant (squared Planck length), 
is the tension (or energy density per unit length) of the string, and units c = h = 1 
are used. In view of the significance of high- dimensional space-times in various aspects 
[241 123] . we consider a generalization of a cosmic string to a codimension-2 brane in 
(d + l)-dimensional space-time with squared length element 

ds 2 = dt 2 - dr 2 - v- 2 r 2 d<p 2 - dz 2 d _ 2 , (8) 

where r and cp are the polar coordinates of the conical surface, 

(j = 3, d) are the Cartesian coordinates of flat (d — 2)-dimensional space, and 

v= (1 -AG^y 1 . (9) 
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The transverse size of the cosmic string is neglected, and both the spatial and the 
bundle curvatures are nonvanishing and singular (as distributions) in the codimension-2 
brane (i.e. point in the d = 2 case, line in the d = 3 case, surface in the d = 4 case and 
(d— 2)-hypersurface in the d > 4 case). Thus, the cosmic string background in which the 
scalar matter field is quantized is characterized by two dimensionless parameters: AGfi 
and e$(27r) _1 , where e is the coupling constant of the scalar to the gauge field forming 
the string and should not be mistaken for the electromagnetic coupling constant. 
The d'Alembertian in the cosmic string background takes form 

n = d?-A^-J2d], (10) 



3=3 



where 



A( tr ) = r^drrdr + v 2 r~ 2 (d v - i^j (11) 

is the two dimensional (transverse) Laplacian and dj = ^j. Green's function in the 
cosmic string background obeys equations 

(D x + m 2 )G (t, r, <p, z d _ 2 ; if, r' , <p', z' d _ 2 ) 
= (D x > + m 2 )G (t, r, <p, z d _ 2 ; t' , r', <p', z^_ 2 ) 

= 5(t - t') V -5{r - r')A(<p - V ')5^ (z d _ 2 - z' d _ 2 ) , (12) 

where A(<p) = (27r) _1 e%nip i s the delta-function for the angular variable, A(<y? + 27r) = 

A((/?), Z is the set of integer numbers. It should be noted that, apart from the overall 
phase factor, Green's function is periodic in the value of gauge flux $ with the period 
equal to 27re -1 , i.e. it depends on 

Z7T Z7T 

where [it] denotes the integer part of quantity u (i.e. the integer which is less than or 
equal to u). Imposing physically plausible conditions on Green's function at small and 
large distances from the string, we solve (12) unambiguously and obtain the following 

\ The limit of the transverse size tending to zero may cause problems which were for a long discussed 
in the literature SHI EH [SD] . It was claimed that the source (stress-energy tensor) of space-time 
(8) makes sense as a distribution [!)TJ [55]. However, a class of regular metrics can provide for the 
distributional source which is concentrated on a codimension-1, and not higher, brane [47] . Various 
limiting procedures and wider classes of metrics were considered for the case of a codimension-2 brane, 
yielding an ambiguity in the relation between the energy density per unit length of the string and 
the deficit angle of space-time [13 HE1 HH]- Finally, a class of metrics was found providing for the 
distributional source with support on a codimension-2 brane and with the ambiguity eliminated [50] . 
This class comprises metrics which are either vanishing (in the case of the negative deficit angle 
2tt(1 — v~ r ) < 0) or diverging (in the case of the positive deficit angle, < 2n(l — v^ 1 ) < 2ir) 
along a normal to the brane in its vicinity. Thus, parameter v in (8) can take values in the range 
< v < oo. 
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expression (see appendix): 

G (t, r, ip, z d _ 2 ] t', r', <p', z d _ 2 ) = 



im (d-l)/2 



(27r)( d+1 )/ 2 



cxp 



ie$ 
~2n 



(d-l)/4 



^(d-i)/2 - t') 2 + (r - r') 2 + (z d _ 2 - z^_ 2 ) 2 + 4rr'sin 2 \fa - <p' + 27r/)(2z/)-^ 

{-(t - t') 2 + (r - r') 2 + (z d _ 2 - z^ 2 ) 2 + 4rr'sin 2 [(y9 -</?' + 2tt/)(2z/)- 1 ] } 

oo 

~ J fif (u, up - iff) 



X ■ 



K( d -i)/2 mJ-(t - t') 2 + (r - r') 2 + ( Z(i _ 2 - z' d _ 2 ) + 4rr' cosh 2 («/2) 



-(t - V) 2 + (r - r') 2 + (z d _ 2 - z^_ 2 ) 2 + 4rr'cosh 2 (M/2) 



(d-l)/4 



(14) 



where the summation is over integers / that satisfy condition 

V?' - g _ v_ j V' -V v_ 
2tt 2 2tt + 2' 



(15) 



K a (u) is the Macdonald function of order a, and 



1 [ e i ( 1 - i? )^+^) cosh(Fz/w) - e -^+^) C osh[(l - F)vu) 



cosh(zz-u) — cosh((/9 + un) 
e i(i-F)^-un) cos h(Fz/w) - e -^(v-^) cosh[(l - F)Au] 1 
cosh(zz-u) — cos((/9 — vk) J 

It should be emphasized that the coincidence limit, x' — > x, is to be understood as 
limit tp' —> tp taken before the coincidence limit of at least another one variable, say, 
before r' — > r. In this case the singularity of (14) (term corresponding to I — 0) at 
x' — > x coincides with that of Green's function in the Minkowski space-time 

im (d-l)/2 



(17) 



G (t, r, <p, z d _ 2 ; r', <p', ^ d - 2 j\^ ~ ( 2n y d+ i)/2 
tf(d-i)/ 2 (my/-^ - t') 2 + (r - r') 2 + (z d _ 2 - z^ 2 ) 2 + 4rr'sin 2 [(y9 - uj>)/2\) 



x 



{-(t - t') 2 + (r - r') 2 + (z d „ 2 - z^ 2 ) 2 + 4rr'sin 2 [(v9 - ip')/2}} 



(d-l)/4 



note that the leading singularity is independent of mass: 

'd-l\ 1 



4 7r (d+l)/2 J 



x — X' 



(here T(-u) is the Euler gamma function). Thus, omitting the term corresponding to 
I — in (14) in the coincidence limit corresponds to subtracting (17) in the coincidence 
limit from (14), and this renormalization procedure yields finite (renormalized) vacuum 
expectation values. 



Induced vacuum energy-momentum tensor 



4. Cosmic string background: renormalized vacuum energy-momentum 
tensor 

The / 7^ terms in the finite sum in (14), as well as the last term given by an integral in 
(14), contribute to the renormalized vacuum expectation value of the energy-momentum 
tensor. Its calculation is simplified by noting the following properties: 

(d t + d t ,)G(x; x') = 0, 

{dj + d r )G(x\ x') = 0, (18) 
{d v + d^)G(x; x') = 0. 

In the coincidence limit we get 

lim d t G(x; x')t R ) = lim djG(x; x')tm = 0, (19) 

x'-tx x'-¥X 

and, consequently, the following off-diagonal components of the induced vacuum energy- 
momentum tensor are vanishing: 

Due to the next relation in the coincidence limit, 

\im(d r -d r/ )G(x;x') {R) = 0, (21) 

the remaining off-diagonal component is vanishing as well: 

t 12 = 0. (22) 

Further, since the Klein-Gordon equation is satisfied by Green's function at x ^ x', 
the following relations are valid: 



(<fc,V£V£ - m 2 ) G(x; x') = - 



d r d r i + ^ 1 d r rd r + r' 1 d r >r'd r ^ 



V^G(x; x') = g w \^-d 2 + J2 d j + r ~ 1,9 rr<9 r - m 2 j G(x; x'). 

In view of the above, we get 

D 00 (x; x') = -d 2 - 2 - ^ (d r d r , + r- l d r rd r ) , 

D u (x; x') = \d r d r , - \d 2 r + 2 (t - ±) r-'dr, 



D 22 (x; x') = v 2 r 



2-2 



d t + J2 d J + 2 + \) r'^rrdr + 2 - d r d 7 



D jj (x; x') = -d 2 + 2 - ^ {drdr' + r~ l d r rd r ) . 



G(x; x'), (23) 
(24) 

(25) 
(26) 

- - iir , (27) 
(28) 
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Applying (25)- (28) to the regular part of Green's function according to (7) and taking the 
coincidence limit, we obtain the renormalized vacuum expectation value of the energy- 
momentum tensor: 

+ (29) 

where 

4 /mx (d+i)/2^ / /vrV 3 ^)/ 2 riV ln\- 2 

Q°o(r) = 5}(r) = -p^ ( 7 ) E cos ( 2 ^) ( sin 7 J { 



+2(1-40 



K(d+i)/2 (zmr sin ^ - 2(1 - 4£)mr sin l -^-K( d+3 )/ 2 ^2mr sin ^ | , 



(30) 



rf« = - psp (7) E cos(2/F7r) ( sm 7) 



x 



. /7T . 

sin — - 4£ 



-2 



^(d+i)/2 ( 2mr sin ^ ) , (31) 



4 , TO x(*-i)/2!"? / /vrV 3 ^ 72 

^ = -(SjW I7J 4. cos ( 2 ^) (an ~ ) 



. In 

sin — - 4£ 

v 



1=1 

K{d+\)/2 2mr sin — - 2mr sin — A( rf+3 )/2 2mr sin — 



v V ^ / 

16z/ /mx^+i)/ 2 7 , „(3-d)/2 



,(32) 



00 



s ° (r) = s ^ r) = (47)7^ (7) J dvA{v] v) W=i 

1 

x {[v- 2 + 2(1 - ^)]K {d+1)/2 (2mrv) - 2(1 - 40mr^X (d+3 )/ 2 (2mr^)} , (33) 

00 

*\<r) = (47r) (, + 3)/2 (7) J dvA (v, F, v)-j==(y- 2 - ^)K {d+m (2mrv), (34) 

1 

00 

s ' (r) = (47)(^)7 (7) J dvA{v > V) W=1 {V - 40 

1 

x [K {d+1 y 2 (2mrv) - 2mrvK {d+3)/2 (2mrv)] , (35) 
and 

A(v; F, v) = -fi£(2arccosh?;, 0) 

sin(Fz/7r) cosh[2(l — F)z/arccosht>] + sin[(l — F)vir] cosh(2Fz/arccoshf ) 

cosh(2z/ arccoshf ) — cos(u it) 

(36) 

clearly, tensor q£ is zero for < v < 2. 
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One can verify that relation 

d r t\{r)+r' 1 [t{(r)-t 2 2 (r)] =0 
is valid; hence, the induced vacuum energy-momentum tensor is conserved: 

*A = °- 

Taking the trace of the tensor, we get 

I"/2] 



9 

(37) 
(38) 



(47r)( d+1 )/ 2 \ rJ 



—mr sin — K 

v 



m 




In 


— 


sin 




r 




V 


(d+ 


3)/2 


(> 



m\ v-f , nln N / /vr^ (1_c ° /2 
— j 2^ cos(2/Fvr) I sin — 

i=i <J 

In 



In 



K {d+i)/2 ( 2mrsin — 



2mr sin 



?7T 



V 



- m K(4-i)/ 2 I 2mrsin 



(39) 



32u 



(4vr)(rf+3)/2 V r 



mx (d-l)/2 /• ^ !,(l-«0/2 

dt>A(t>; F, i/ 



= {m c -o 



m 



x [ir (d+1)/2 (2mrf ) - mrvK (d+3)/2 (2mrv)] - m 2 K {d _ 1)/2 {2mrv)} , (40) 
where 

H c = {d-l){Ad)-\ (41) 

Hence, the trace is proportional to m 2 at £ = £ c , and conformal invariance is achieved 
in the massless limit at £ = £ c : 



lim g^D* 



m— s-0 



0. 



(42) 



that is in accordance with a general treatment in [44j |45] . 
5. Various limiting cases 

Using the asymptotics of the Macdonald function at large values of its argument, 
we obtain the asymptotics of the induced vacuum energy-momentum tensor at large 
distances from the cosmic string. Namely, we get in the case < v < 2: 

t\(r) = -(2mr)~ l t 2 2 (r), mr > 1, (43) 

where 

2u 



*8(r) = t]{r) = t 2 2 (r 



n , cos[(F-l/2)i/7r] 

( 47r )(d+i)/ 2 ^ V sin ( z/vr /2) vW 



x exp(— 2mr){l + 0[(mr) ]}, mr > 1 
in the case < v < 2, and 

t° (r) = t](r) = t 2 2 (r) 



(44) 



2 , , /m\ d / 2 

(4n)^ {1 -^ m U 



cos(2Fvr) + (2F - 1 



sin(2F 



n 



y/nmr 



x exp(— 2mr){l + 0[(mr) ]}, mr > 1 



(45) 
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in the case v — 2. We get in the case 2 < v < oo: 

ttx -i 



t{(r) 



2mr sin — 
v 



mr 3> 1, 



(46) 



cos(2F7r) (sin — 



7T\ -d/a 



7T\ 2 

1 - 4£ ( sin - 



m\ d / 2 
m [ — ) 

r / 



[I -40m (-) 



" (47r)rf /2 

x exp (— 2mr sin — ) {1 + O^mr)^ 1 ]}, mr > 1, (47) 

*8(0=<J(r) = ^^008(2^) (sin-) 

x exp (— 2mr sin —J {1 + 0[(mr) -1 ]}, mr > 1. (48) 

Asymptotics (44) in the case d = 3 was obtained earlier in [37] . It should be noted that 
the asymptotics in the case 2 < v < oo, (46)-(48), is less decreasing than that in other 
cases. The asymptotics at small distances from the cosmic string coincides with the case 
of the vanishing mass of the scalar field; this case will be considered in the next section. 
In the case of a global string ($ = 0), we get 

H2] , ,_ x (3-d)/2 



t" (r)\ F= o = t>(r)\ F=0 



+2(1-40 



K 



(d+l)/2 



(47r)(d+i)/2 V r / 

?7T 



m\ (d+l)/2 



£( sin v) 



Z=l 



sm 



2mr sin 



v 



Itt 



- 2(1 - 4£)mr sin —K (d+3)/2 



2mr sin • 



8z/ sin(z/7r) /m\( d + 1 )/ 2 



(47r)(d+3)/2 V r y 

v (3-d)/2 



dv 



x 



cosh 2 (z/ arccoshu) — cos 2 (z/7r/2) 
2(1 - 40mruir (d+3) /2(2mru)} , 



{[tr 2 + 2(l-40]iW)/ 2 (2mrt;) 



(49) 



*}(r)| 



F=0 



( 47r )(d+i)/2 V r y 



sm ■ 



v 



4£ 



Mx ,u, [f/2] / , N (3-d)/2 

z=i ^ v ' 



In 



K( d +i)/2 2mrsin 



+ 



8i/sinfi/7r) /m\(<i+i)/2 



(47r)(^+3)/2 



di> 



v^^l 



,(3-d)/2( v -2_ 4 ^ 



cosh (i/arccoshf) — cos 2 (z/7r/2) 



%+i)/2(2mra), 



(50) 



*2(r)| 



F=0 



(4 vr )(d+i)/2 V r y 



m x (d+l)/2 



(3-d)/2 
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-2 



X 



sin — - 4£ 



v 



lir\ Itt / /7r 

^V(rf+i)/2 ( 2mr sin y I - 2mr sin — A(d +3 )/2 I 2mr sin — 



^ 8z/sin(i/7r) /mx^ 1 )/ 2 ^ dv 



( 47r )(d+3)/2 V r 7 y v^^T 

i 

v (3-d)/2f v -2 _ 4 ^ 

X co S hVarcco S H-co B VV2) [ A '<«" 2 < 2m ™> " ^^(^(a™™)] ■ («) 

In the case of the vanishing string tension (/i = 0) when the string is reduced to a 
merely magnetic one, the result of [40J is recovered: 

oo 

n/ ,-, 16sin(Fvr) /m\(<*+i)/2 f dt> . 

«r)U = «r)U = p^( 7 ) y - 7 _cosh[(2F-l)arccosht;] 

i 

{[ v -2 + 2(1 _ 40]X (d+ i)/ 2 (2mr ? ;) - 2(1 - 4£)mrvK {d+s)/2 (2mrv)} , (52) 

oo 

16sin(Fvr) /m\(^+i)/2 f dw , r/n „ 1N 

*i(r)U = ( 7 ) y cosh[(2F - l)arccosh,] 

i 

xv (l-J)/2( v -2 _ A ^ K(d+l)/2 (2mrv), (53) 

oo 

a , l6sin(Fvr) /m\W /" dw , r/n „ 1N , . 

^Wl-i = (47r) (d + 3)/2 J y "y^T cosh [( 2F - l)arccosh^] 

i 

xv (i-«0/2( v -2 _ 4£) [if (d+ . 1)/a (2mrv) - 2m™# (( ^ )/2 (2m™)] ; (54) 

it should be noted that temporal component (energy density) (52) at £ = 1/4 was 
obtained in I4T1I421. 



6. Massless quantized matter 

Using the asymptotics of the Macdonald function at small values of the argument, we 
obtain the induced vacuum energy-momentum tensor in the case of the massless scalar 
field: 

^' lm =° = n^&m r ' d ' 1 {[P(i(F ' v) ~ ^Pd-2{F, v)\ dia g (i, i, -d, i, . . . , i) 

+4(£ - QP d - 2 (F, z/)diag(rf - 1, -1, d, d - 1, . . . , d - 1)} , (55) 
where £ c is given by (41), and 

P d (F, u) = P^(F, v) + P { d s) (F, v), (56) 

, , v-d-i 

Pi q) (F, = cos(2/Fvr) sin - , (57) 

i=i \ u J 

oo 



1 
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The tensor trace is determined by Pd-2{F, v) only: 

9^'\ m=0 = iEC^ r-^ - Qd(d - l)P d _ 2 (F, i/). (59) 

In the case of the vanishing string tension, the integral in (58) is taken, yielding 
(see 00]) 

W sin(F 7 r)2^r(^ + F)r(^ + l-F) 

Pd (F ' 1} " — T(dT2) ' (60) 

while Pj 9) (F, 1) = 0. 

In the case of the nonvanishing string tension, the integral in (58) is explicitly taken 
for odd values of d only [36]. For this task one should use the integral representation 
through contour C= on the complex w-plane, see figure 2 and, e.g., (A. 18) in appendix, 

P?(F, = d ^W ^M^. (61) 

87Ti J suih(ra /2J 

c= 

If d is odd, then the last factor in (61) gives a pole at w = 0, and contour C= is deformed 
to encircle this pole, as well as other poles on the imaginary axis; hence the integral 
takes form 

if>(F, ,) = -P«(F, „) + -L I^B^F)^ (oddrf), 



n £ i 
n>0 



(62) 

where the origin is circumvented anticlockwise, and B n (u) is the Bernoulli polynomial 
of order n (see, e.g., [53]). Only the simple pole contributes, yielding 

(d+l)/2 

P d {F,u)= C d + l v 2n B 2n (F) (oddd), (63) 

n=0 

where 

C'n 1 = 77t~tt A /dw^-^-sinh 2 ^^)]-^ 1 )/ 2 . (64) 



(2n)! 4vri 

The integral in (64) is in a familiar way transformed into an integral over contour C = 
resulting in 

CI+ 1 = - J- J- / d WW 2 "- 1 [-sinh 2 («;/2)]-( d+1 )/ 2 
2n (2n)! 4vri 7 L v ' n 



1 1 



c= 

t J du [(u + ivr) 2 "" 1 -{u- i7r) 2n ~ 1 ][cosh(u/2)]- d - 1 . 



(2n)! 2yri 
o 

Although the integral in the last line of (65) can be taken with the help of 



(65) 



d\ \ (2n)\ \ 2 

k 



i ^ (-i)^ 2 v (-mk+dy. { 



n{2n-l)j^ (2/)! k\[k + (d + l)/2] 2 ("-0-i 



k>() 
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the obtained expression is hardly operative for n > 1. The more efficient expressions 
can be obtained for separate consecutively decreasing values of n: 

2 d 



(jd+l _ ^_-Q(d+l)/2 
C d+1 = (_-Q(d-l)/2 



2 d - 2 d + 1 



(d- 


1)! 


6 




-4 


d+1 


(d- 


3)! 


72 


2 d- 


-6 


d+1 


(d- 


5)! 


432 



C£} = (- 1 ) M/2 7T L ^T ^ ( j + rf J > (67) 

-/--. < ') .. , -„ ,,, _ 35 + 5 + 3 

and so on; also, expressions in terms of the generalized Bernoulli numbers can be used, 
see [36]. In view of the above, we get the following results for Pd{F, v) (63) up to, for 
instance, d — 9: 

P 1 (F, v) = \- v 2 B 2 (F), (68) 
6 

Ps(F, u) = ± + jB 4 (F) + ^P X {F, u), (69) 

PAF, u) = + — BAF) - —BAF) + -PAF, u), (70) 
5V ' ; 9450 15 y J 45 y J 5 V ' h K J 

oon Q, .4 A j, 6 ..8 a 

P 7 (F, v) = -=^— + —BAF) - —BAF) + —BAF) + -PAF, u), (71) 
V ; 198450 315 V ; 189 V ; 315 V ; 7 V ; V ; 

12637 19z/ 4 , . 14z/ 6 . . v 8 , . 2z/ 10 . . 

PAF, v) = + BAF) BAF) + BAF) B 10 (F) 

JV ' ; 22453200 3240 V ; 1215 V ; 405 V ; 14175 V ; 

+Ip 7 (F, u), (72) 
where the explicit form of relevant Bernoulli's polynomials is 
B 2 (F) = l-F(l- F), B A {F) = -1 + F\l - Ff, 

B e (F) = 1 - l -F\l - Ff - F\l - Ff, 

Bs(F) =~ + \f\1 ~ F) 2 + V(l - Ff + F\l - Ff, (73) 

B W (F) = 1 - ^F 2 (l - F) 2 - 3F 3 (1 - Ff - V(l - Ff - F 5 (l - Ff. 

Note the appearance of the quadrupled conformal coupling as a factor in the last terms 
in (69)-(72) (i.e. 4£ c at d = 3, 5, 7 and 9, respectively); due to this circumstance, the 
tensor in the conformally invariant case (£ = £ c ) does not contain terms which are 
linear in v 2 and in F(l — F), only the higher integer powers (from z/ 4 to v d+l and from 
[F(l — F)] 2 to [F(l — F)] (d+1) / 2 ) occur. This fact holds true for arbitrary odd values of 
d, being a consequence of relation 

dC d+1 = (d-l)C d -\ (74) 
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which follows from the explicit expression given by (66) at n = 1 

r d+l _ r (~2~) / 7 r^ 

Ca """fifny (75) 

As to coefficient C$ +1 , it is positive and decreasing from 1/6 to as d increases from 1 
to oo, 

oo 

Co +1 = / -^- 2 [cosh(u/2)]- d -\ (76) 
J It -|- 7T 


It should be noted that relations (68) and (69) are known for quite a long time, see 
[351136]. 

If d is even, then the integrand in (61) has a branch point singularity at w = and 
representation (58) is more convenient for a numerical evaluation of Pff (F, v) in this 
case. 



7. Summary 

In the present paper we have shown that a cosmic string induces the finite energy- 
momentum tensor in the vacuum of the quantized massive scalar field. Although a 
cosmic string in astrophysical applications is characterized by small positive values of 
the deficit angle of conical space, that corresponds to the range of parameter v (9) 
restricted to < 1 - u' 1 < 4 x 10" 7 [9], our consideration extends to the most general 
cases of conical space with all possible positive and negative values of the deficit angle, 
that corresponds to < v < oo (see footnote f). Another generalization which is 
employed in the present paper is that of a cosmic string to a codimension-2 brane in the 
static (ultrastatic) locally flat (d + l)-dimensional space-time, see (8). The transverse 
size (thickness) of the brane is neglected, whereas the mass of the quantized scalar field 
is taken into account and the coupling of the field to the scalar curvature of space-time 
(£) is assumed to be arbitrary. 

The general expression for components of the induced vacuum energy-momentum 
tensor is given by (29)-(36). A remarkable feature is a possibility of analytic continuation 
to complex values of space dimension, that yields the tensor components as holomorphic 
functions of d on the whole complex d-plane. In the case < v < 2, all the dependence 
on the brane parameters (tension and flux) is contained in function A (36), whereas the 
dependence on d is factored otherwise. The Aharonov-Bohm effect [55] is manifested 
in the dependence of (30)-(35) on the fractional part of e$(27T) _1 (13) rather than 
on the flux itself. Moreover, the tensor components are symmetric under substitution 
F — > 1 — F; consequently, the induced vacuum energy-momentum tensor is even under 
charge conjugation. 

A completely novel feature is that the tensor components in the case 2 < v < oo 
attain an additional to (33)-(35) contribution which is given by (30)-(32). This feature 
is analogous to the finding of the authors of [56j [571 EE] that the incoming wave in the 
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problem of quantum-mechanical scattering on a conical singularity is a superposition of 
a finite number of distorted plane waves propagating in different (rotated) directions. 
To be more precise, the number of the additional terms in the tensor components, that 
equals {u/2\ (see (30)-(32)), is related to the number of waves propagating out of the 
shadow or double- image region in conical space, that equals 2[z//2] + 1 (see [59, 60J). 

Finally, it should be noted that the tensor components decrease exponentially at 
large distances from the brane, see (43)- (48). If the mass of the quantized scalar field is 
zero, then expressions for the tensor components are simplified considerably, see (55)- 
(58). In the case of odd dimension of space, the dependence on the flux is contained 
in Bernoulli's polynomials of even order, which are multiplied by v in the power of the 
same order; we give explicitly the results for the cases d — 3, 5, 7 and 9, see (68)-(72). 
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Appendix 

We start with the kernel of the resolvent of operator — A^ tr ^ + m 2 : 

R"(r, cp; r', cp') = (r, ^\(- ^ + m 2 - u 2 y l \r\ cp'), (A.l) 

where u is a complex parameter. The resolvent kernel in the cosmic string background 
obeys equations 

(-A^+m 2 -Lo 2 )R u (x; x') = (-Aj, r) +m 2 -w 2 )/T(x; x') = -6(r-r')A(<p-<p'),(A.2) 

where x = (r, ip), x' = (r',ip') and A^ is given by (11). We impose the 
condition of regularity of (A.l) at small distances, r — > or r' — > 0, while the 
condition at large distances is that (A.l) behave asymptotically as an outgoing wave: 
r -1 / 2 exp(ir\/uj 2 — m 2 ) at r — > oo or r'~ 1 / 2 exp^r'y/cu 2 — m 2 ) at r' — >■ oo, where a 
physical sheet for square root is chosen as < ArgV^ 2 — m 2 < n (Inii/u 2 — m 2 > 0). 
As a result we get 

R W (r, ip; r', <p') = — ^ { d ( r ~ r ')^[r i -e$(27r)-i](^)4|n-e<I»(27r)-i|(^ / ) 

+9(r' - r)/ !/ | n _ e$(27r) - 1 |(«r)X I/[n _ e#(27r) -i ] («r / )} e in ^'\ (A.3) 

where I a (u) and K a (u) are the modified Bessel functions of order a (see, e.g. [53]), 6{u) 
is the Heavyside (step) function, and 

k = -Woo 2 - m 2 , -tt/2 < ArgK < tt/2 (Re k > 0). (A.4) 
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C- 



71 



w 



C + 



-71 



Figure 1. The vertical part of contour C + is infinitesimally close to the imaginary axis 
from the right; the vertical part of contour C_ is infinitesimally close to the imaginary 
axis from the left. 



Choosing r' > r, we rewrite (A. 3) as 



v 



i? w (r, ip; r , <p ) = — exp 



Y J I ^-F){Kr)K v{k _ F) {Kr')^ k ^ 



fe>i 



+ ^2lv(k>+F)(Kr)K Hk , +F) (K,r')e lk ' 



k'ez 

k>>0 
oo 

f 

x / — exp 

J y 



where 



— exp 



9 / 9 i /9 

1^ lyi" _l_ ty' ^ 



2y 



2 , /2 \ 



fcez 
fe>i 



fc'ez 

fc'>0 



(A.5) 



(A.6) 



F is the fractional part of e$(27r) 1 , see (13). 

Using the Schlafli contour integral representation for the modified Bessel function 

I a ( y ) = — [ dwe ycoshw ~ aw = — — / dwe ycoshw+aw , 
2ni J 2ni J 

c+ C- 

where contours C + and C_ in the complex w-plane are depicted in figure 1, we perform 
summation in (A.6) to get 



v 



dw c ycosh w ~ Fuw |- ([ w e v coshw ~ Fuw 



]_ g— vw— up 



C- 
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i 


t | w 

— — 


i 




■J 


i 


• : 


>- 

D 




N — "* — 



Figure 2. Two straight horizontal lines and circles enclosing simple poles at w = 
— i(<£ + 2irl)i>~ 1 {— 7r <Imii) < 7r). 



- / 

2717 J 
C++C- 



dw e v coshw - FL,w . 



1 



Q — VW — lip 



(A.7) 



Further, the union of contours C + and C_ can be continuously deformed into a contour 
depicted in figure 2. Taking the residues of poles on the imaginary axis, we get 

y cos [(<£> + 2nl)u x ] } 



2ni 



dw e ycoshw 



e (l-F)(uw+i<p) 



(A.8) 



(A.9) 



where the summation is over integers / satisfying condition 

(p V if v 

and contour C= consists of two straight horizontal lines depicted in figure 2. Substituting 
(A.8) into the last line of (A. 5), we get 



R?(r,<p; r',<p') = — exp 



x 



e 2ilFn K (k^J (r - r') 2 + 4rr'sm 2 [(</? - + 2nl)(2v)- 1 ]) 
i ^ 

'')] / I \ 

— K I Ky (r — r') 2 — 4rr'sinh (w/2) I 



v 

27Ti 



e {l-F)[vw+i(<p-<p')\ 
1 — e vw+i(<p-<p r 



(A.10) 



here the summation is over integers / satisfying condition (15). 

With the use of the resolvent kernel one can obtain Euclidean Green's function in 
2+1-dimensional space-time with the Wick- rotated time axis, t = — ir: 

G E (r, r, w r', r', y?') = (r, r, (-«9 2 - A (tr) + m 2 )" 1 |r', r', 
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(A.11) 



Substituting (A. 10) with u = ip into the last line of (A. 11) and performing the 
integration over p, we get 



G (r, r, <p; t', r', <p') = exp 



2tt 2 



ie$ ,' 



x 



E e * 



2ilFn 



exp |-m^/(r - r') 2 + (r - r') 2 + 4rr'sin 2 \(p - p' + 27r/)(2z/T 1 ] } 



/ du> 

2vri 



v/(r - t') 2 + (r - r') 2 + 4rr' sin 2 [(y? - ^ + 2tt/)(2z/)- 1 ] 

e (l-F)[(^+i(^-^')] eXP 



-m-y/ (t — t') 2 + (r — r') 2 — 4rr'sinh 2 (w/2) 



c= 



(r - r') 2 + (r - r') 2 - 4rr'sinh 2 (w/2) 



. (A.12) 



Euclidean Green's function in 3+1-dimensional space-time with imaginary time 
t = — ir can be obtained as a double Fourier transform of the resolvent kernel: 

G E (r, r, p, z; t', r', p', z>) = (r, r, p>, z\ (-<9 2 - A (tr) - d 2 z + m 2 )" 1 |r', r', p' , z') 



^ J dqe iq{z - z,) J dpe ip(r - T ' ) J R i v / « 5 +^(r, ip; r', p'), 



(27T) 

— oo 

or equivalently: 



(A.13) 



oo 

^ / 



G*(r, r, ^ z; r', r', p', z') = — / dqe^-^G{r. r. p: r ,'. r 

Substituting (A.12) with m — > a/ m 2 + q 2 into (A. 14), we get 

~ie$ 



>^/m 2 +q 2 



(A.14) 



G K (r, r, </?, z; r', r', <//, ^ 



m 

W) 



exp 



2tt 



(p - p') 



x 



K x (my/(r - t') 2 + (r - r') 2 + (z - z') 2 + 4rr'sin 2 [(v? - p> + 2nl)(2u)- 1 ] 



v 

27ri 



dw 



y/(r - r') 2 + (r - r') 2 + (i - z') 2 + 4rr' sin 2 [{p - <p' + 2tt/)(2z/)- 1 ] 



X — e vw+i(ip-ip') 



C= 



\J{j - t') + (r - r') 2 + (z - z') 2 - 4rr'sinh 2 (w/2) 



(A.15) 



Proceeding in this way one can get Green's function in space-time of arbitrary 
dimension, since 

G e (t, r, p, z 1} z d _ 2 ; t', r', p', z[, . . . , z' d _ 2 ) 



= 27 / ^ q(Zd - 2 ' Z ' d - 2) G E {r, r, <p, z d _ 3 ; r', /, <p', z[, . . . , ^j^^— 
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(2vra)( rf - 2 )/ 2 



i / Tfi 

%-2)/2(aV? 2 + ™ 2 ) = 
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(A.16) 

m \(d-i)/2 K {d _ 1)/2 (mVa 2 + z 2 ) 



( a 2 + z 2)(d-l)/4 > 

(A.17) 



which is derived with the help of the integral representation for the Macdonald function 



(see, e.g., [53] 



oo 

K a (2bc) = \ Q° J dyy^exp - c 2 ^ . 

o 

Thus, Euclidean Green's function in d + 1-dimensional space-time is given by expression 

m (d-l)/2 



G E (r, r, <f, z d _ 2 ; r', r', ip', z^_ 2 ) 



■ cxp 



ie$ 



,2i/F7r 



(d-l)/2 



t — r') 2 + (r — r') 2 



( 27r )(«M-l)/2 

+ (z d _ 2 - z d _ 2 ) 2 + 4rr'sin 2 [(^ - <p> + 2tt/)(2z/)- 1 ]) 



{(r - r') 2 + (r - r') 2 + (z d _ 2 - z' d _ 2 ) 2 + 4rrsin 2 [(<^ - <// + 2vr/)(2z/)- 1 ] } 



in (d-i)/4 



^ r e (l-F)[vw+i(<p-ip')] 

2ti\ J ^ W 1 - e ^ +i ^-^') 



A' 



(d-l)/2 



X 



) 2 + (r — r') 2 + (z^_ 2 — z^_ 2 ) 2 — 4rr'sinh 2 (u>/2) 



m\ [t — T 



[(r - r') 2 + (r - r') 2 + (z d _ 2 - z d _ 2 ) 2 - 4rr'sinh 2 (w/2)] 



(d-l)/4 



(A.18) 



Transforming the integral over contour C= into an integral over the real positive 
semiaxis, we get 



G E (t, r, y>, z d „ 2 ; r', r', yj', z^_ 2 ) 



m 



(rf-l)/2 



(27r)( d+1 )/ 2 



exp 



ie$ 



2i/F7r 



X 



(d-i)/a {m\/(r- r') 2 + (r - r 



'\2 



+ (z d _ 2 - z^_ 2 ) 2 + 4rr'sin 2 [(^ - ^ + 2tt/)(2z/)- 1 ]) 



x 



{(r - r') 2 + (r - r') 2 + (z d _ 2 - z^_ 2 ) 2 + 4rr'sin 2 [(^ - <p' + 2vr/)(2z/)- 1 ]} (d 1)/4 

oo 


A (d _ 1)/2 (my/ (r - r') 2 + (r - r') 2 + (z d _ 2 - z^_ 2 ) 2 + 4rr'cosh 2 (u/2) 



[(r - r') 2 + (r - r') 2 + (z d _ 2 - z^_ 2 ) 2 + 4rr'cosh 2 (w/2)] 



(d-l)/4 



(A.19) 



where VL F {u, <p) is given by (16). Going over to real time, we obtain Green's function 
(14). 
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